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ABSTRACT 

We propose that when neutron stars in low-mass X-ray binaries accrete sufficient mass 
and become miUisecond pulsars, the interiors of these stars may undergo phase transi- 
tions, which excite stellar radial oscillations. We show that the radial oscillations will be 
mainly damped by gravitational- wave radiation instead of internal viscosity. The gravi- 
tational waves can be detected by the advanced Laser Interferometer Gravitational- Wave 
Observatory at a rate of about three events per year. 

Subject headings: dense matter — gravitation — stars : neutron 
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1. INTRODUCTION 



Gravitational wave astronomy may soon become an observational science, since three 
gravitational wave experiments, including the Laser Interferometer Gravitational- Wave 
Observatory (LIGO) (Abramovici et al. 1992), are under construct. In astrophysics, neu- 
tron stars are widely believed to be the most promising source of gravitational radiation 
(for detailed reviews see Thorne 1987, 1995), which from collapse of the cores of massive 
stars may provide a signature for the features of supernova explosion (Burrows & Hayes 
1996), which from starquakes of pulsars may reveal the physics of the stellar interiors 
(Zimmermann & Szedenits 1979), and which from mergers of binary neutron stars may 
give information about the equation of state (EOS) of nuclear matter at high densities 
(Shibata, Nakamura & Oohara 1992, 1993; Rasio & Shapiro 1994; Zhuge, Centrella & 
McMillan 1994; Davies et al. 1994; Ruffert, Janka & Schafer 1996). In this Letter we pro- 
pose a new possible origin of gravitational- wave bursts. We argue that when neutron stars 
in low-mass X-ray binaries accrete sufficient mass and become millisecond pulsars the in- 
teriors of these stars can undergo phase transitions, which excite stellar radial oscillations, 
producing strong gravitational wave bursts. 



2. THE MODEL 



2.1. Evolution of Neutron Stars in Low-Mass X-Ray Binaries 



According to the standard scenario of evolution of low-mass X-ray binaries (Bhat- 
tacharya & van den Hcuvel 1991), mass is transferred from the companion to the neutron 
star, which is spun up to a millisecond period. The mass-transfer from the companion 
keeps driving the two stars apart while processes such as orbital gravitational radiation or 
magnetic braking keeps driving the two stars closer. These processes keep the system in a 
steady mass-transfer state throughout the evolutionary timescale of the companion while 
the accretion rate of the neutron star is near the Eddington value. Thus, the neutron 
star can accrete mass > 0.5Mq in ~ 10^ years and become a millisecond pulsar (van den 
Heuvel & Bitzaraki 1995a, b). If we make an assumption that the masses of neutron 
stars before accretion are 1.4M0, which is supported by the current theories of Type II 
supernova explosion and observations of masses of pulsars (e.g., the Hulse-Taylor binary 



2 



system), then the stars in an evolutionary timcscalc > lO^yrs must become rather mas- 
sive ones (> I.8M0). Now we ask a question: what could possibly occur in the interiors 
of these massive neutron stars? 

To study this question, we first analyze possible EOSs for neutron stars. So far there 
have been many approaches to determine an EOS for dense matter through the many- 
body theory of interacting hadrons. Unfortunately, these approaches have given EOSs 
with different stiffnesses and in turn very different structures of neutron stars. However, 
the EOSs should be constrained by the observations as follows. First, Link, Epstein & 
van Riper (1992) used a model-independent approach to analyze the post glitch recovery 
in four isolated pulsars (Crab, Vela, PSR 0355+54 and PSR 0525+21) which are likely 
to be isolated IAMq neutron stars, and concluded that soft EOSs at high densities are 
ruled out. More detailed analyses of the post glitch curves of the Crab and Vela pulsars 
also draw similar conclusions (Alpar et al. 1993, 1994). Second, if the EOSs in cores 
of neutron stars with mass ~ 1.4M0 were soft, the massive compact objects after the 
accretion phase of low-mass X-ray binaries could be black holes (Brown 1988). In fact, 
these objects have been identified as millisecond pulsars. This means that soft EOSs are 
unlikely to occur in neutron stars with mass ~ 1.4M0. For these two reasons, we can 
assume that the EOS in neutron stars with mass ~ IAMq is moderately stiff to stiff. 

The above assumption is consistent with recent theoretical studies of the EOS for 
dense matter at high densities. First, because of the strong repulsion between nucleons 
and nucleon holes in the spin-isospin interaction, pion condensation is unlikely in neutron 
stars (Brown et al. 1988; Baym 1991). Second, the possibility of kaon condensation in 
dense matter was suggested by Kaplan and Nelson (1986), whose basic idea is that the 
energy of a negative kaon is lowered by interaction with nucleons. In neutron star matter 
in beta equilibrium one expects negative kaons to be present if the energy to create one 
kaon in the matter is less than the electron chemical potential. However, Pandharipande, 
Pethick and Thorsson (1995) studied kaon- nucleon and nucleon-nuclcon correlations in 
kaon condensation in dense matter and found the kaon energy is much larger than that of 
the previous studies (sec, e.g.. Brown et al. 1992; Brown et al. 1994; Thorsson, Prakash 
& Lattimer 1994; Lee et al. 1995; Thorsson & Wirzba 1995). On the other hand, if we 
adopt electron chemical potentials of the calculations of the simple parametrized models 
of Prakash, Ainsworth and Lattimer (1988) which represent a number of more realistic 
models and compare these potentials with the kaon energies of the Hartree calculations 
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(Pandharipandc et al. 1995) in Figure 1, then we can see that the density for kaon 
condensation is ~ (5-7)po (where po is the nuclear density), which significantly exceeds 
the central densities of neutron stars with mass ~ IAMq and with moderately stiff or stiff 
EOSs. Several authors (e.g., EUis, Knorren & Prakash 1995; Schaffner & Mishustin 1996; 
Dai & Cheng 1996) have used the relativistic mean-field approach to study kaon energies 
in neutron star matter, and have also drawn similar conclusions. Third, it is thought that 
the density for deconfinement of nuclear matter to two-fiavor quark matter is near (6-9) po 
(Baym 1991). Therefore, the recent theoretical studies of dense matter also indicate that 
pion (or kaon) condensation or quark matter is unlikely to occur in neutron stars with 
mass ~ IAMq. 

We now turn to the study of what happens in neutron stars with moderately stiff 
to stiff EOSs when they accrete mass > 0.5Mq, viz., when their masses increase from 
IAMq to more than I.SMq. It is clear that the central densities of these massive stars 
are ~(5-7)po- Once this condition is reached, several physical processes will take place 
in the stars. When e~ — > K~ + z/g, kaon condensation occurs spontaneously in the 
stellar interior. The appearance of this new phase destroys the hydrostatic equilibrium 
of the star due to the softening of the EOS of the core matter. This implies a structural 
transition of the whole star into a new, stable configuration with a condensed core of the 
new phase. The structural transition — neutron star corequake — occurs on the dynamic 
timescale of milliseconds, in which the interior temperature may increase to ~ lOMeV 
both due to the rapid compression of dense matter (Haensel, Denissov &: Popov 1990) 
and due to the energy released by reaction e~ — > K~ + v^- Alternatively, if the central 
nuclear matter is deconfined into two-flavor quark matter, the quark matter will convert 
to three-flavor quark matter because strange matter is thought to be more stable than 
nuclear matter. Thus a strange matter seed is formed in the interior, and subsequently 
the strange matter will begin to swallow the neutron matter in the surroundings. This 
process should proceed in a timescale ~ tens of milliseconds due to a detonation mode, 
and the interior temperature increases to > 10 MeV because the chemical energy of the 
two-flavor quark matter is dissipated into thermal energy (Dai, Peng & Lu 1995). The 
conversion of neutron stars to strange stars has been suggested as a possible origin of 
cosmological 7-ray bursts (Cheng & Dai 1996). Furthermore, the phase transitions of 
massive neutron stars to either stars with kaon condensation cores or strange stars can 
stimulate stellar radial oscillations. In next subsection we will show that these newly born 
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compact objects may be a strong source of gravitational-wave radiation if their rotation 
periods are of the order of miUiseconds. 



2.2. Damping of Radial Oscillations 



Radial oscillations are damped not only due to dissipation of the vibration energy of 
stellar matter into heat but also due to conversion of this energy into gravitational-wave 
radiation. We first focus on the case in which a massive neutron star undergoes a phase 



neutrinos are trapped in the interior and form an ideal Fermi-Dirac gas with /^^ ^ kT 
(where T is the interior temperature) because the neutrino mean free path is much less 
than the stellar radius for KT ~ 10 MeV. Subsequent evolution of the newborn star, which 
is analogous to that of a protoneutron star formed from supernova explosion (Burrows 
& Lattimer 1986), can be divided into three stages. We neglect gravitational radiation. 
The first stage is deleptonization, whose timescale {rai) is of the order of Is (Sawyer & 
Soni 1979), in which the interior temperature may increase significantly. In this stage the 
radial oscillations are damped through the following reaction 

p + e~ ^ n + Ug, (1) 

where p and n represent quasiprotons and quasineutrons respectively. Using the reaction 
matrix element of Brown et al. (1988), we have derived the net rate per baryon for this 
reaction at nonequilibrium. 



where Gp is the weak coupfing constant, 9c is the Cabibbo angle, 6 is the chiral angle 
for kaon condensation {9 < 60°) (Brown et al. 1994), qa is the Gamow-Teller coupfing 
constant, m* and m* are the effective nucleon masses, and ^ = /ip + /le — /J'n — IJ'u 
with /li being the chemical potential of particle i. In deriving equation (2), we have 
assumed kT/c <^ <^ |P^ with P,^ and being the neutrino and neutron momenta 
respectively. As the steps shown by Dai & Lu (1996), using this reaction rate, we have 
further derived the bulk viscosity as 



transition to a star containing a kaon condensation core. After e 



K + z/g, electron 



1 



cos^ 9c cos'' {9 /2){1 + 3^^ )^*^*^^^2^[^2 ^ ^27TkTf] , (2) 
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771 ~ 1.8 X 10^^cos-2(^/2)i; 



IMeV 



g cm ^ s 



(3) 
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where Y^, Yy and Yn are the particle concentrations, and p is the stellar density. Therefore, 
the damping timescale (Sawyer 1980) is given by 

Tvl = ^PR^^I' (4) 

where R is the stellar radius. Inserting equation (3) into this equation, we obtain 
c 0.52 cos^(^/2)FV3F-/3r-/3^2 j j \^ 

where is in units of 10^ cm. For the typical particle concentrations and p ~ 6po, 
R ~ 10^ cm and kT - lOMeV, r„i - 12rV3y-2/3y^-4/3 ^ > i2s > r^i. 

In the second stage, the stellar interior has practically no trapped lepton-number ex- 
cess, as compared with catalyzed matter. The diffusion of v^V^ is then driven by the 
temperature gradient. Locally, the equilibrium distribution function of v^V^ can be ap- 
proximated by the Fermi-Dirac one with a zero chemical potential. The neutrino diffusion 
timescale {Td2) is of the order of 40 s (Sawyer & Soni 1979). Using the analogy of deriving 
the bulk viscosity of (Haensel & Zdunik 1992), we obtain the damping timescale 

r,2 ^ 1.8 X 10^ cos\e/2)Y-^'^Rl f j s , (6) 
Clearly, r„2 > Td2- 

Third, after the neutrino diffusion, the temperature decreases to ~ 1 MeV, and neutri- 
nos can escape freely from the star. The damping timescale becomes 

- - 2.6 X 10-cos-W2)n-"»i;^ [ff {^)' s, (7) 

where uj is the oscillation frequency. Therefore, we can conclude that the radial oscillations 
are not damped very efficiently by the bulk viscosity until the stage at which neutrinos 
escape freely. 

We now consider gravitational radiation from rapidly spinning and oscillating neutron 
stars. The timescale for this process (Chau 1967) is 

r,^0.41M,:X-'(^)' s, (8) 
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where Mi.g is the stellar mass in units of 1.8Mq, and P is the stellar rotation period. Here 
we have assumed that nucleons in the stellar interior are nonrelativistic and degenerate, 
and thus the adiabatic index is equal to 5/3. When Tg < Tdi + Td2, viz., 

P<6.2f]/XWms, (9) 

where = {rai + Td2)/4:0s, the gravitational radiation can damp the radial oscillations 
very efficiently. The frequency of the gravitational waves is euqal to u; = 27r/r, where r is 
close to 5 X 10~^ s for a typical neutron-star model (Glass & Lindblom 1983) (corrected for 
the gravitational redshift) with mass M ~ I.SMq and radius R ~ 10^ cm. The strength 
of the waves can be estimated using the quadrupole approximation to the Einstein field 
equations (Thorne 1987). This approximation shows that the gravitational strain is given 

by 

^ ~ 1.5 X 10-3^^ ^ P ^ 



(o.l) (2ms) (lOOMpc) ' 



c4 r " \6poJ VO.iy V2msy V^OOMpc 

where Q is the source's quadrupole moment (Chau 1967), r is the distance of the source 
from Earth, and a is the relative oscillation amplitude (= 5R/R). Particularly, in the 
case of phase transitions of neutron stars with stiff EOSs to stars with kaon condensation 
cores, a is expected to be about 0.1. Furthermore, the characteristic gravitational strain 
is 

. .vs . 4.3 X nri^i (4) (^) (^) (^) (H) 

where n is the number of cycles of gravitational waves in the duration ~ Tg. The observed 
number of low- mass X-ray binaries in our Galaxy is ~ 10^ (van Paradijs 1995), and thus 
the rate for phase transitions of their neutron stars is ~ 10~^yr~^ because the typical 
accretion timescale is ~ 10^ years. An alternative estimation based on the number of 
millisecond pulsars and their hfe time also gives a similar value (Gheng & Dai 1996). 
Therefore, the rate detected by the advanced LIGO detector (Abramovici et al. 1992) is 
estimated to be about three events per year. 

Alternatively, neutron stars in low-mass X-ray binaries may accrete sufficient mass to 
convert into strange stars, as suggested by Gheng & Dai (1996). Next we would discuss 
this case. First, the timescale for damping radial oscillations due to bulk viscosity is > 10 s 
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for a high temperature > lOMeV (Madsen 1992; Dai & Lu 1996), so the gravitational 
radiation damping mechanism is also more efficient in this case. Second, for the conversion 
of neutron stars with stiff EOSs to strange stars, the relative oscillation amplitude {a) 
may not be less than 0.1. Third, the density for deconfinement may be larger than that 
for kaon condensation, and thus the oscillation frequency of strange stars is larger. If 
gravitational waves discussed in this work are observed by the advanced LIGO in the 
future, then observations, in principle, can distinguish between these two kinds of phase 
transitions. In addition, gravitational radiation for the case in which neutron stars convert 
to strange stars is hkely to occur together with cosmological 7-ray bursts, because the 
fireballs formed during the conversion have very low baryon contamination (Cheng & Dai 
1996). 

3. CONCLUSIONS 

We in this Letter have suggested that during the evolution of neutron stars in low- 
mass X-ray binaries the stars may undergo phase transitions to stars containing kaon 
condensation cores or strange stars when these neutron stars accrete sufficient mass from 
their companions. The phase transitions can excite stellar radial oscillations, which pro- 
duce strong gravitational wave bursts if the stellar rotation periods are of the order of 
milliseconds. The study of such gravitational radiation may provide information for the 
high-density EOS and the physics of phase transitions. In addition, the rate detected by 
the advanced LIGO is estimated to be about three events per year. 

Gravitational waves from colliding neutron stars must have a continuous spectrum, 
but waves from phase transitions of neutron stars appear to have a delta-function like 
spectrum. This signature is expected to be confirmed by very near future observations . 

We thank W.M. Suen for bringing this problem to our attention, W.Y. Chau for dis- 
cussions on gravitational radiation, and Samuel Wong and M.C. Chu for discussions on 
equations of state at high densities. K.S.C. thanks a RGC grant of Hong Kong for support; 
Z.G.D. thanks the National Natural Science Foundation of China for support. 
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FIGURE CAPTION 



FIG. 1. Energy of a single negative kaon in neutron star matter and electron chemical 
potential as functions of density. The three solid lines are kaon energies from the Hartree 
calculations of Pandharipandc ct al. (1995) for square wells of radii R = 1 and 0.7fm, and 
for a Yukawa potential. The dashed lines labeled by PALI, PAL2 and PAL3 represent 
electron chemical potentials calculated from the parametrized models of Prakash et al. 
(1988) corresponding to three different forms of F{u) which parameterizes the potential 
contribution to the symmetry energy. 
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